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I The most precise current determination of the Unitarity Triangle angle (p^ is presented. 

fSJ ' It is obtained using the decay B'^ D^*^K'^ , with Dalitz plot analysis of the subsequent 

Zj'. D ^ KsTT^TT- decay. The result is 03 = 68° l|^o(stat) ± 13°(syst) ± 11° (model). 

1 Introduction 

The year 2004 has seen a great number of important new results from Belle. Studies of time- 
^ ! dependent asymmetries in i? ^ tt+tt^ revealed the second observation of CP violation, and the 
^ I first evidence for direct CP violation, in the B system ^J. More recently, further evidence for 
^ direct CP violation has been found in i? — *^ K~7r~^ decays; this effect is also seen by BaBar 
Q ■ The decay B —>■ 7r°7r° is now observed with more than 5a [3], and evidence for B p^n^ has 
. been seen jl]. Studies of new particles observed at Belle have stimulated theoretical interest [S], 
0-) ! and further excitement has been aroused by new results on the time-dependent asymmetries in 
JlLiI b ^ s penguin transitions jE]; new results on 5 — > KsKgKs and B — > Ksir^'j were shown for 
^ ' the first time in FPCP2004 [?]. To complete this brief selection of the many other new results 
from Belle jHj, the first studies of the forward-backward asymmetry in B ^ K*l^l~ have been 
performed [Oj. 

This plethora of new results is made possible by the excellent performance of the KEKB 
\ accelerator. KEKB is an asymmetric-energy e^e~ (3.5 GeV on 8 GeV) collider ^Hj, operating at 
the T(4S') resonance {^/s = 10.58 GeV) with a peak luminosity that exceeds 1.3 x 10^^ cm~^s~^. 
At the time of FPCP2004, the total integrated luminosity was approaching 300 fb~^ (at the 
time of writing it is approaching 330 fb~^), however the results presented here use the data 
sample of 253 fb~^ accumulated by summer 2004, corresponding to approximately 275 x 10^ 
BB pairs. 

The Belle detector is a large-solid-angle magnetic spectrometer that consists of a three- 
layer silicon vertex detector (SVD), a 50-layer central drift chamber (CDC), an array of aerogel 
threshold Cerenkov counters (ACC), a barrel-like arrangement of time-of-flight scintillation 
counters (TOF), and an electromagnetic calorimeter comprised of CsI(Tl) crystals (ECL) lo- 
cated inside a superconducting solenoid coil that provides a 1.5 T magnetic field. An iron 
flux-return located outside of the coil is instrumented to detect mesons and to identify 
muons (KLM). The detector is described in detail elsewhere ^llj. Two different inner detec- 
tor configurations were used. For the first sample of 152 million BB pairs, a 2.0 cm radius 
beampipe and a 3-layer silicon vertex detector were used; for the latter 123 million BB pairs, 
a 1.5 cm radius beampipe, a 4-layer silicon detector and a small-cell inner drift chamber were 
used 
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2 Principles of 03 Extraction from D^*>K 

The main objective of the B factories is to make precise measurements of the elements of the 
Cabibbo-Kobayashi-Maskawa (CKM) quark mixing matrix [13j, As is well known, the unitarity 
of the CKM matrix results in the expression 

Vu,v:, + v,,v:, + Vuv:, = o, (i) 

which can be represented as a triangle in the complex plane (the so-called Unitarity Triangle). 
Two popular (and various other less popular) sets of names for the angles of this triangle are 
used, 



01^/5 



are 



a 



are 



7 = arg 



(2) 



and 



It has long been known that decays of the type B DK are sensitive to 03 
various methods using different final states of the neutral D meson have been proposed 
These are based on two key observations: neutral and mesons can decay to a common 
final state, and the decay B^ D^*'>K^ can produce neutral D mesons of both flavours via 
h — > cus and h ucs transitions, with a relative phase 6'_|_ that depends on both strong and 
weak phase differences between the two amplitudes, ^+ = 5 + 03. For the charge conjugate 
mode, the relative phase ifi 6_ = 5 — 03. The size of possible CP violation effects depends on 
the ratio of the magnitudes of the interfering amplitudes, which we denote by r. 

Recently, three body final states common to and D , such as i^'^vr+Tr" ^E], were 
suggested as promising modes for the extraction of 03. The principle of the technique is 
as follows (for more details see [^). We define the amplitude for Kst^^t^~ at a 

point on the Dalitz plot given by (m 



to be f{m\,rn?_)-, where m\ and m?_ are the 



squared invariant masses of Ks'ir'^ and Kstt , respectively. Assuming CP invariance in D me- 

,2 ^2 ^ rjnj^g amplitude for 



son decay, the corresponding amplitude for decay is then /(m^,mj 



B~^ — >■ DK^ can then be written as /(m 



-re^'^^^'f'^^ f{m^_,m\), while that for B~ 



DK- 



is f{m'^_,m\) + re*'''^~"'^^^/(m^, m^). Thus, once the D decay model / is fixed, the values of r, 
03 and 5 can be determined by simultaneously fitting data from B^ and B~ decays. 

The same arguments also apply to similar decays such as B^ D*K^ . Some care must be 
taken, since for each B decay mode the values of r and 5 can be different. It has recently been 
pointed out that for B^ —>■ D*K~^ there is an effective shift of vr in the value of 6 between the 
cases that the decays D* Dtt^ and D* — ^ D'j are used 18j. 

The D decay distribution can be measured at B factories using the large samples of flavour 
tagged D mesons produced in —>■ D^tt^ (and charge conjugate) decays following continuum 



e^e 



cc reactions. However only f{m 



m^) can be measured resulting in some uncertainty 
due to the model assumed to describe the phase variation across the Dalitz plot. 

This technique was first implemented by Belle, using a 140 fb~^ data sample, to obtain 
03 = 77° ll^o(stat) ± 13°(syst) ± 11° (model) HH). Recently the same method has been used by 
BaBar, with a data sample of 211 x 10*^ BB pairs, to obtain 03 = 88° ± 41°(stat) ± 19°(syst) ± 
11° (model) 1201. 

The results presented below are preliminary updates of the previous Belle analysis. More 
details can be found in 
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3 Dalitz analysis 



3.1 Extraction of D Decay Model 

We select candidates for the D — > Ksn^n^ decay, and combine with slow charged pions (vr^) 
to make D*^ candidates. We require D*^ candidates to have momenta above 2.7 GeV/c^ in 
the centre of mass (cm) frame to reject background from B decay. We fit the mass difference 
AM = Mp.^^+^^^± — Mk gTT+TT- to obtain signal and background fractions, and then require 
144.6 MeV/c^ < AM < 146.4 MeV/c^. The distribution of background events across the 
Dalitz plot is obtained from a sideband region of AM. We then fit the candidates to a model 
including the components listed in Table. ^ with efficiency, resolution and background taken 
into account. The result of the fit is shown in Fig. ^ The /ndi of the fit is 2543/1106, 
however fine tuning the model has little effect on 03 and possible discrepancies are taken into 
account in the model uncertainty. 
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Figure 1: From left to right: Dalitz plot distribution of tagged D decays; projections onto the 
, ml and _ variables. 
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3.2 Selection of B Decay Candidates 

We reconstruct neutral D mesons in the D —>■ Ks^'^tt' mode, and D* mesons by making 
combinations with 7r° candidates. We then combine these with prompt tracks, which are se- 
lected as kaon candidates on the basis of information from the CDC, TOF and ACC systems. 
The selection of B candidates is based on the CM energy difference AE = J^Ei — -Ebeam and 

the beam-constrained B meson mass Mbc = \/ E^cam ~ i^PiY^ where -Ebeam is the CM beam 
energy, and Ei and Pi are the CM energies and momenta of the B candidate decay prod- 
ucts. The requirements for signal candidates are 5.272 GeV/c^ < Mbc < 5.288 GeV/c^ and 
\AE\ < 0.022 GeV. The signal and background fractions are obtained from binned fits to the 
AE distributions, which are shown in Fig. |21 The signal is modelled with a Gaussian, an addi- 
tional Gaussian is used to describe the background from 5^ D'^*'>tt^ with the prompt pion 
misidentified as a kaon, that peaks at AE ~ 50 MeV. The remaining background is modelled 
by a linear function. For B^ — ^ DK^ the selection efficiency is 11%, there are 276 candidates 
satisfying all selection criteria of which 209 ±16 are found to be signal events. For i?^ D*K^ 
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K*(1410)-7r+ 


0.35 ±0.04 


253 ±6 


0.14% 


Ko*(1430)+7r- 


2.15 ±0.04 


348.7 ± 1.1 


7.4% 


ir*(1430)-7r+ 


0.52 ±0.04 


89 ±4 


0.43% 


A'2*(1430)+7r- 


1.11 ±0.03 


320.5 ± 1.8 


2.2% 


K*(1430)-7r+ 


0.23 ± 0.02 


263 ±7 


0.09% 


K*(1680)+7r- 


2.34 ± 0.26 


110 ±5 


0.36% 


A'*(1680)-7r+ 


1.3± 0.2 


87± 11 


0.11% 


nonresonant 


3.8 ±0.3 


157 ±4 


9.7% 



Table 1: Results of the fit to determine the D decay model. 

the selection efficiency is 6.2%, there are 69 candidates satisfying all selection criteria of which 
58 ± 8 are found to be signal events. The Dalitz plot distributions of candidate events are also 
shown in Fig. |21 The Dalitz plot distributions of the different background sources are obtained 
using sideband regions in data or Monte Carlo simulation. 

3.3 Extraction of 03 

In order to test the analysis procedure, we use control samples of D^*'>'k^] these have 

similar topology and phenomenology to our signal, but in this case the ratio of amplitudes is 
expected to be small (r ~ 0.01). We obtain results which are consistent with expectation, and 
take possible discrepancies into account in the systematic error. 

We then perform unbinned maximum likelihood fits to the D^*'>K^ samples. In order 

to extract the three free parameters 03, r and 5, we fit the B^ and B^ samples simultaneously. 
However, to display the result, as in Fig.|21 we show likelihood contours in (Re(re*^±), Im(re*^±)), 
that are obtained separately for B^ and B~ . CP violation is seen as a difference between the 
B^ and B~ contours. 

Since the value of r is positive definite, and the uncertainty on ^3 depends on r, we do not 
rely on the errors obtained from the likelihood curves |22] 1 but instead use a frequentist approach 
to obtain confidence regions, which are also shown in Fig. |21 The results for B^ —>■ DK^ are 

03 = 64° ± 19°(stat) ± 13°(syst) ± ll°(model), 
r = 0.21 ± 0.08(stat) ± 0.03(syst) ± 0.04(model), (3) 
5 = 157°± 19°(stat) ± ll°(syst) ±2r(model). 
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Figure 2: (Top) AE distributions and fit results for (left) DK^ and (right) D*K^ 

candidates. (Bottom) Dalitz plots for (from left to right) B^ DK^, B^ DK^, B^ 
D*K^ and B~ —>■ D*K~ candidates. The axes are exchanged between B^ and B^ candidates 
so that CP violation effects appear as differences between the B~^ and B~ distributions. 



which correspond to r > at 99.3% CL and CP violation at 94% CL. The results for B^ 
D*K^ are 

03 = 75°±57°(stat)±ll°(syst)±ir(model), 
r = 0.12l[]:J?(stat) ±0.02(syst) ±0.04(model), (4) 
5 = 321° ± 57°(stat) ± ll°(syst) ± 21°(model), 

which correspond to r > at 56% CL and CP violation at 38% CL. The dominant systematic 
uncertainty is due to possible fit bias (note again that more details can be found in ^^). 

To combine the results from both modes we use a frequentist approach with Feldman- 
Cousins ordering and obtain 

03 = 68°l}^o(stat) ± 13°(syst) ± ll°(model), (5) 

which corresponds to CP violation at 98% CL. The two standard deviation interval including 
systematic and model uncertainties is 22° < 03 < 113°. 



5 




Figure 3: Likelihood contours for (left) — > DK^ and (right) — > D*K^ candidates. In 
each pair, the leftmost plot is the result of the unbinned fit, and the rightmost is the outcome 
of the frequentist treatment, as described in the text. 



4 Summary 

The most precise measurement of 03 is presented. The value is consistent with the previous 
measurements from Belle jTHj and BaBar [20] • We obtain a significantly smaller statistical error 
than BaBar due to the larger value of r obtained (there is, however, no significant discrepancy 
in the values of r obtained from these and other measurements ^^). 
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